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Abstract 

By using Jet Calculus as a consistent framework to describe multiparton dynamics 
we explain the peculiar evolution equation of fracture functions by means of the recently 
introduced extended fracture functions. 
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1 Introduction 



Fracture functions KJ have been introduced to interpret within the framework of perturba- 
tive QCD semi-inclusive deep inelastic processes in the target fragmentation region. The 
formulation of these processes by means of fracture functions does allow to extend the inter- 
pretation in terms of QCD-improved parton model to the complementary region of target 
fragmentation. This implies a description of deep inelastic processes in terms of the novel 
fracture functions in addition to the usual factorizable structure and fragmentation functions 
convoluted with hard point-like cross sections. 

The introduction of these new objects requires however the formulation of an additional 
factorization hypothesis whose validity has been recently argued on the basis of the cut 
vertex formalism together with the use of infrared power counting. These results have been 
obtained within the theoretical framework of (0 3 )6 held theory [0, [| and have been later 
confirmed in QCD in Ref. {§. 

In order to show that the new factorization hypothesis is correct, in Ref.[[J new objects 
have been defined called extended fracture functions, which depend also on the momentum 
transfer t between the incoming and outgoing hadron. By an explicit one-loop calculation it 
has been verified that these object do factorize and it has been also observed that they show 
a new logarithmic dependence on the ratio of the scales Q 2 and £ . As a result extended 
fracture functions obey an evolution equation different from the one followed by ordinary 
fracture functions. 

The aim of this note is to show in a simple and direct way that extended and ordinary 
fracture functions are closely connected. As a consequence the corresponding evolution 
equation do follow one from the other. We work within the framework of Jet Calculus by 
applying the corresponding rules to the evolution equations. The set of Jet Calculus rules 
allows the interpretation in terms of QCD-improved parton model of the results obtained 
by using a different approach ||. By using this method we can give in fact a definition of 
the extended fracture function in the region where t is a perturbative scale ( g ^ 1). As 
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a consequence a DGLAP evolution equation for extended fracture functions follows and the 
inhomogeneous evolution equation for the ordinary fracture function is recovered. 



2 Evolution pattern 



In a deep inelastic semi-inclusive reaction, a hadron A with momentum p is struck by a far 
off-shell spacelike photon with momentum q and a hadron A' with momentum p' is inclusively 
observed in the final state. Let us define as usual 

, = g (!) 

and choose a frame in which p = (p + ,p_, 0) with p + p_ and pq ~ p+q~- 

As far as we keep away from the target fragmentation region the cross section factorizes 
as follows 

*j = J ^r^F^x'^Q^a^x/x' \z/Z \Q 2 )D A ,(~z> \Q 2 ) (2) 

where F A (x, Q 2 ) and D A ,(z,Q 2 ) are the structure and fragmentation function respectively, 
&ij(x, z, Q 2 ) is the hard semi-inclusive cross section and we have defined z = pp' /pq ~ p'_/q~ 



Eq. (D) expresses the fact that, as long as the produced hadron has a large transverse 
momentum p'^ ~ Q 2 (i.e. z is finite) it can be thought as a product of the fragmentation of 
the active parton. However in the last few years, especially after the appearance of diffractive 
deep inelastic events at HERA, particular attention has been payed to hadron production in 
the target fragmentation region (i.e. the region z — > 0), where eq. (Q) fails. In Ref. [U it has 
been proposed for target fragmentation the additional factorized term 

a T = j ^M; A ,(x',z,Q 2 )a,(x/x',g 2 ) (3) 

where z = p'q/pq ~ p' + /p+ represent the momentum fraction of the hadron A' with respect 
to A. Here M AA ,(x, z, Q 2 ) is the fracture function, giving the probability of finding a parton 



i with momentum fraction x in the hadron A while another hadron A' with momentum 
fraction z is detected. 

The fracture function is expected to satisfy an inhomogeneous evolution equation [0, 
and this fact has been verified at one-loop level in Ref . || . 

Let us consider the case in which the momentum transfer t = — (p — p') 2 <C Q 2 is also 
measured. The current and target fragmentation contributions are in this case 

oj= / -j-tF 1 a (x', Q 2 ) aij(x/x', z/x'z r , tx'/z', Q 2 ) D A ,(z', Q 2 ) (4) 

J X z 

and 

r dr' 

a T = J — M AA ,(x',z,t,Q 2 )a t (x/x',Q 2 ) (5) 

The function Ai A A ,(x, z,t,Q 2 ) is the extended fracture function. In Ref. it has been 
shown that this object can be defined just in terms of a new cut vertex. Whereas M\ A , (x, z, Q 2 ) 
is expected to satisfy an inhomogeneous evolution equation |TJ], the extended fracture func- 
tion obeys a simple DGLAP evolution equation 0. 

This facts can be understood by making the following observations. In the region Aq CD <C 
t <C Q 2 the hard semi-inclusive cross section in eq. ([D develops large logQ 2 /t corrections 
which have to be resummed ||. Instead of absorbing these logs in we choose to move 
them in the extended fracture function which, by using Jet Calculus, can be defined in the 
perturbative region of t as 

M AA ,(x,z,t,Q 2 ) = ^ J F A (w,t)dwP t kl (u)duD lA ,(^t)d£ 
x E J k (r, t, Q 2 ) dr5(x — ruw) 5(z — w(l — u)£) 

= a^f d_.^h_ pA pu JL f z \ ■ 2 
Ztti J rww — x/r rw \w — x/r J 

•1-2 



2nt 



("-( ^—,F A (w, t)Pt l (-) D lA i t) E{{x/r, t, Q 2 ) 

(6) 
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where the integration limits have been obtained implementing momentum conservation. 
We work within the leading logarithmic approximation. Here P-{u) and P/ fc (u) are regular- 
ized and real Altarelli-Parisi vertices || respectively, and the function E{(x,Q 2 ,Q 2 ) is the 
evolution kernel from Qq to Q 2 , which obeys the DGLAP evolution equation: 



., , -P>(u)E (x/u,QlQ<). 



(7) 



Eq. (^) is represented in Fig. 1. The evolution kernel (the black blob in Fig. 1) resums the 




Figure 1: Extended fracture function in the perturbative region oft 



dependence on Q 2 and t in terms of logarithms of the form log ^- and at first order in 



has the following expression 

El(x,t,Q 2 ) ~ 6}6{l-z) + g^'(z) log 2. 



(8) 



In this way we have absorbed these large logarithms in the definition of Ai l A A ,(x, z, t, Q 2 ) 
which in the perturbative region of t is resolved in the convolution of the objects contained 
in the dotted box in Fig 1. 
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From eq. ([|) it appears that the Q 2 dependence of the extended fracture function is 
completely given by the evolution kernel E^(x, Qq,Q 2 ) and therefore the evolution equation 
is 

Q 2 -^M\ A ,(x,z,t,Q 2 ) = ^p- ^-Pi(u)M\ A ,(x/u,z^Q 2 ) (9) 



1-z 



at least in the perturbative region of t, i.e. t ^ ^-qcd- Therefore the anomalous dimen- 
sion associated to the extended fracture function is the same that controls the evolution of 
ordinary structure functions. 

Let us suppose now that, as argued in Ref. @, the evolution equation @ applies also 
within the small t region and define the integrated fracture function as an integral over t up 
to a Q 2 -dependent cut-off of order Q 2 , say eQ 2 , with e < 1 



(10) 



M j AA ,( Xj z,Q 2 )= [ dtM\ AI (x,z,t,Q 2 ). 
Jo 

By taking the logarithmic derivative of eq. fllCf ) we have that 

Q 2 ^M{Ax, z, Q 2 ) = ^ Q dtQ 2 ^M 3 AA ,(x, z, t, Q 2 ) + eQ 2 M AjA ,(x, z,eQ 2 , Q 2 ) 

' y * {Q ~ ] f Q dt f —Pi(u)M AA ,(x/u,z,t,Q 2 )+eQ 2 M AA ,(x,z,eQ 2 ,Q 2 ) 
Jo J^- u 

1 — 2 

f -Pttu)MX A ,(x/u,z,Q 2 )+eQ 2 M AA ,(x,z,eQ 2 ,Q 2 ). (11) 

Z7T / x U 

We see appearing an inhomogeneous term in the evolution equation of M° A A ,(x, z, Q 2 ) which 
arises from the Q 2 dependence of the upper integration limit. This inhomogeneous term 
depends on the value of M° A A ,(x, z, t, Q 2 ) in the perturbative region of t where eq. (||) 
holds. By using the boundary condition 

Ei(x,Q 2 ,Q 2 ) = 515(1 -x) (12) 



2tt 

a s (Q 2 ) f 1 du 
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on eq. (P) we get for the inhomogeneous term, up to log e corrections, 



eQ 2 M' AA ,(x,z,eQ 2 ,Q 2 ) = 

2tt A r Jz+r w \ w ~ r ) \w — r J 

27r J z+x w(w — x) \W/ \W — X J 

«S(Q 2 ) f~ z d,U U , „2\A?// Nr, Z' ZU /ton 



By substituting eq. (|T3| ) in eq. ([TT]) the evolution equation for z, Q 2 ) appears 

to be 

Q 2 ^Mi A ,(x,z,Q 2 ) = ^Pi(u)MX A ,(x/u,z,Q 2 ) 
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Z7T Jj. X[l — U) \X{1 — U) J 

that is exactly the one proposed in Ref.0. 

Furthermore the perturbative definition ([]) can be used to derive the evolution with t 
of the extended fracture function. In order to make the notation less cumbersome, it is 
convenient to define 

1 r-l-z 

dzz m dxx n M J AA ,(x,z,Q 2 ) = M^i A \Q 2 ). (15) 
Jo 

One can verify that eq. (^|) becomes 

MUt, Q 2 ) = ^P^Fi +n (t)D k m (t)Ei\t, Q 2 ) (16) 
where we have defined || 

P^= f duu m {l-uYPl\u) (17) 

J 

and 



El l (t,Q 2 )= [ duu n Ej(u,t,Q 2 ) 
Jo 
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The t evolution equation contains several terms: the first comes from the canonical scale 
dependence of the extended fracture function, the second from the scale dependence of 
as, and the third from the scale dependence of the evolution kernel. There are then two 
inhomogeneous terms which follow from the t dependence of structure and fragmentation 
functions respectively. By explicitly deriving eq. (|16D with respect to t we get 



t^Mi n (t,Q 2 ) = - (6™(l + p Q a s (t)) + ^-Af) M p mn (t,Q 2 ) 

— i— a s(^ pkU MP pp (t) D k (t) E jl (t O 2 ) 
(27r) 2 t m m+n n 

+ ^^PlK +n (t)A k j;D? m (t)El\t,Q 2 ) (19) 

where 

Ag= Cduu n Pi{u). (20) 



We conclude this section by stressing that one can define the fracture function in a more 
general way as an integral up to an arbitrary Q 2 -dependent integration limit 

rt 2 (Q 2 ) 

M AA ,(x,z,Q 2 ) = / dtM> AA ,(x,z,t,Q 2 ) (21) 
Jo 

In this case the evolution equation reads 



Q 2 -^M AAI (x,z,Q 2 ) = -Pl[u)MX AI {xlu,z,Q 2 



l-z 



2vr t 2 (Q 2 ) J x r J r r(l - u) 
x D ltA . (-r^s, UQ 2 )) E{ (x/r, t 2 (Q 2 ), Q 2 ) (22) 



or, by taking double moments, 

Q 2 ^- 2 MUQ 2 ) = ^^^MUQ 2 ) 



(23) 



In the case in which ^(Q 2 ) = tQ 2 the evolution equation fll4]) is of course recovered. 



One could even define a more general fracture function introducing also a lower integra- 
tion limit ti(Q 2 ) 

ft2{Q 2 ) 

M AA ,(x,z,Q 2 ) = / dtM 3 AA ,{x,z,t,Q 2 ). (24) 

Jt!{Q 2 ) 

and the equation would have an additional inhomogeneous term taking into account the Q 2 
dependence of the lower integration limit 

Q 2 ^ 2 MUQ 2 ) = ^-K p Ml n {Q 2 ) 

+ asit *^ ))Q l^ P™F^ HQ 2 )) Dt [HQ 2 )) El 1 HQ%Q 2 ) 

- as{t ^ ))Q l^ p^F^ hq 2 )) nt HQ 2 )) eS HQ 2 IQ 2 ) ■ 

(25) 

Phenomenologically these two cases would correspond to the production of hadrons 
within the target fragmentation region with transverse momentum below ^(Q 2 ) and be- 
tween ti(Q 2 ) and ^(Q 2 ) respectively. 



3 Sum rules 

In Ref. I] it was shown that fracture functions obey the following momentum sum rule 

I dzzM AA ,(x,z,Q 2 ) = (1 - x)F A (x,Q 2 ) (26) 

A' J 

which accounts for momentum conservation in the s-channel. In this section we want inves- 
tigate if a momentum sum rule holds also in the t-channel. Taking moments of eq. (|14]) with 
m = 0, n = 1 and summing over j we get || 

Q 2 ^J2 M ^ A '(Q 2 ) = ^^F^iQ 2 ) j\uuPt(l-u)D k A '(Q 2 ) (27) 
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where we made use of the well-known property of the splitting function || 

5>? = (2f 



j 

and of the relation 

-i 



o 



du(l-u)y Jf(u) = / du{l-u)Pl(u). (29) 



From eq. (|27| ) it appears that the derivative of 52 . Mq[ (Q 2 ) receives contribution from the 
inhomogeneous term in the evolution equation. This fact suggests that the sum rule could 
be violated because at values of t of order Q 2 the current contribution becomes important 
and the two production mechanisms cannot be disentangled. 

As far as the extended fracture function is concerned, since A4 A A ,(x, z,t,Q 2 ) obeys the 
DGLAP evolution equation, it follows that 

Q2 ~^pYl [ *dxM 3 AA ,(x,z,t,Q 2 ) = (30) 



o 



and momentum conservation in the t channel is recovered. In the perturbative region of t, 
by using eq. (pi), we find 



J^M^'^Q 2 ) = ^F- A (t) f\uuP*{l-u)D k A \t). 



(31) 



4 Summary 

In this note we showed that a formulation of perturbative processes in the target fragmen- 
tation region can be given by using the formalism of Jet Calculus 0. 

A perturbative evaluation in terms of evolution equation and anomalous dimensions of 
extended fracture function is consistent with the evolution equations proposed in Ref.|] 
and |2| for ordinary and extended fracture functions respectively. We showed that the 
inhomogeneous term in eq. (|TJJ) precisely stems from the integration over momentum transfer 
t. 
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Moreover in this formalism the extension to next-to-leading order seems quite natural 
and suggests that eqs. (|J) and QI4]) keep the same structure in the two-loops approximation, 
analogously to what happens for ordinary structure and fragmentation functions |10|| . 

In the region Aq CD <C t <C Q 2 the hard semi-inclusive cross section develops large 
\ogQ 2 /t corrections which need to be resummed. Our approach is based on the idea of 
absorbing such corrections in the definition of A4 AA ,(x,z,t,Q 2 ), namely \ogQ 2 /t are fully 
contained in the evolution kernel E^(x,t,Q 2 ). These logs, which can be treated with stan- 
dard renormalization group techniques, are new and potentially useful to understand the 
dynamics of hadro-production in the target fragmentation region. These logs could affect 
significantly observables like multiplicities and particle distributions and we believe that a 
phenomenological study of the role of these corrections would be useful. 



As originally proposed in Ref. |IJ] and in Ref.|[ll||, there now seems to be a widespread 
consensus that the Q 2 evolution in the diffractive regime can be described in terms of the 
perturbative QCD formalism |3], [12|, The data taken by HI collaboration [14|, showing 



an evidence of a consistent logarithmic scaling violation in diffractive channels have further 
stressed that such an interpretation is not far from the experimental observation. All this 
is consistent with a description in terms of fracture functions. Our results agree with such 
an expectation, although in the limit z — > 1 our formulae will show the appearance of large 
log(l — z) factors which should be resummed as well. For the treatment of such a singular 
region we believe an approach can be used close to the one followed in the inclusive case. 

As a concluding remark we want to stress that the results presented here completely 
agree with those of Ref. [0], obtained within a quite different framework, namely, that of a 
generalized cut vertex expansion which, as a generalization of Operator Product Expansion, 
appears to be founded on a firmer theoretical standpoint. Once more the equivalence of 
parton-like and OPE-like approaches is verified. 



10 



Acknowledgments 



Many of the arguments covered in this letter originate from discussions with G. Veneziano. 
We also thank G. Veneziano for useful comments, and S. Catani, D. Graudenz and J. Kodaira 
for conversations. One of us (M.G.) would like to thank the Department of Physics of the 
university of Florence for the hospitality at various stages of this work. 

References 

[1] L. Trentadue and G. Veneziano, Phys. Lett. B323 (1994) 201. 

[2] M. Grazzini, L. Trentadue and G. Veneziano, Nucl. Phys. B519 (1998) 394. 

[3] M. Grazzini, Phys. Rev. D57 (1998) 4352. 

[4] J.C. Collins, Phys. Rev. D57 (1998) 3051. 

[5] M. Grazzini, Nucl. Phys. B518 (1998) 303. 

[6] K. Konishi, A. Ukawa and G. Veneziano, Nucl. Phys. B157 (1979) 45. 

[7] G. Altarelli, R.K. Ellis, G. Martinelli and S.Y. Pi, Nucl.Phys. B160 (1979) 301. 

[8] D. Graudenz, Nucl. Phys. B432 (1994) 351. 

[9] L. Trentadue, Proceedings of QCD Euroconference 96 (Montpellier, France July 1996) 
published in Nucl. Phys. Proc.Suppl. 54A (1996) 176; see also D. Graudenz, Fortsch. 
Phys. 45 (1997) 629. 

[10] J. Kalinowski, K. Konishi, P.N. Scharbach and T.R. Taylor, Nucl. Phys. B181 (1981) 
221,253. 

[11] A. Berera and D.E. Soper, Phys. Rev. D53 (1996) 6162. 



11 



[12] Z. Kunszt and W.J. Stirling, talk given at Durham Workshop on HERA Physics: Pro- 
ton, Photon, and Pomeron Structure, Durham, England, 17-23 September 1995, and at 
International Workshop on Deep Inelastic Scattering and Related Phenomena (Rome, 
Italy, April 1996), piep-ph/9609245] Z. Kunszt, talk given at Ringberg Workshop on 



New Trends in HERA Physics (Ringberg Castle, Germany, May 1997) H^jZHHHl ■ 



[13] D. De Florian and R. Sassot, [hep-ph/9710205| , to be published in Proceedings of the 
Madrid low-x Workshop (Miraflores de la Sierra, Spain, June 1997). 

[14] HI Collaboration, Z.Phys. C76 (1997) 613. 



12 



